Charge transport in ballistic multiprobe graphene structures 
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We study the the transport properties of multiterminal ballistic graphene samples, concentrating 
on the conductance matrix, fluctuations and cross-correlations. Far away from Dirac point, the cur- 
rent is carried mostly by propagating modes and the results can be explained with the conventional 
semiclassical picture familiar from ray optics, where electrons propagate along a single direction 
before scattering or reaching the terminals. However, close to the Dirac point the transport is due 
to evanescent modes which do not have to follow a rectilinear path. As we show in this Letter, this 
property of the evanescent modes influences the conductance matrix. However, at best it can be 
observed by measuring the cross correlations in an exchange Hanbury Brown-Twiss experiment. 



Transport in undoped graphene is entirely due to 
evanescent modes and can be seen as a form of electron 
tunneling. In contrast, lightly doped graphene supports 
propagating modes that have a linear dispersion relation, 
similar to the k • P approximation of semiclassical elec- 
tron dynamics. Semiclassical electrons can be described 
in terms of classical trajectories and a picture of ray op- 
tics, whereas it is not as clear how evanescent waves move 
in a sheet of graphene. One way to provide more insight 
to this duality of evanescent and propagating modes is 
the study of cross-conductances and cross-correlations of 
electric current in multiprobe graphene structures. This 
is the aim of our present Letter. 

One of the striking properties of electronic transport 
due to evanescent modes in graphene is the "pseudodiffu- 
sive" behaviour in undoped samples, manifesting itself in 
the transmission distribution of ballistic graphene at the 
Dirac point It turns out that all the different cumu- 
lants of current (fluctuations) through such a graphene 
sample behave in the same way as they would in a dif- 
fusive wire 0. However, the nature of conduction in 
the two cases is quite different: in diffusive wires, the 
conduction electrons are propagating but due to mul- 
tiple scattering the information on the propagation di- 
rection in such systems is (almost) lost. One outcome 
of this momentum isotropization can be observed in an 
exchange Hanbury Brown-Twiss cross-correlation exper- 
iment, where the measurement result depends on coher- 
ent processes that connect all terminals involved in the 
measurement Q. Making such a cross-correlation ex- 
periment in conventional ballistic conductors would yield 
a vanishing result, because in the absence of scattering 
propagating modes can only couple pairs of terminals. 

In ballistic graphene, there are no elastic scatterers, 
but due to the evanescent nature of the charge carriers, 
the momentum direction of the electron waves is not well 
defined. Rather, the waves are spread out and a single 
evanescent wave can couple each of the terminals. As 
a result, the exchange Hanbury Brown-Twiss result will 
be very similar to that in diffusive wires (see Fig. [5]). 
The behavior of the two kinds of modes is schematically 
illustrated in Fig. Q] 
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FIG. 1: (Top) A schematic explanation for the difference be- 
tween ray optics and evanescent optics picture of electron 
propagation. Propagating modes (left) couple only two of 
the leads in the absence of elastic scattering. Which lead the 
electron enters depends on the angle of the electron beam at 
the first lead-graphene contact. This fact can be used to esti- 
mate the dependence of the cross-conductances on the spatial 
dimensions of the scattering region (see Eq. (0)). Evanescent 
modes (right) in contrast couple all of the leads. (Bottom) 
Two alternative geometries for the four-probe setup consid- 
ered in this Letter. In setup (a) interfaces between graphene 
and leads 1 and 2 are of zigzag type and interfaces to leads 3 
and 4 are of armchair type. In (b) graphene sheet has been 
rotated by 90 degrees, so that the interface types are switched. 



To compute the conductances and correlators, we em- 
ploy the tight-binding approach, which has been success- 
fully used to describe the transport in graphene @, [E @| 
and is also the starting point in the derivation of the 
Dirac equation describing charge carriers in graphene 0] • 
All the numerical results of this Letter are obtained with 
this approach. On top of this, our aim is to describe 
how well these results can be understood on a qualita- 
tive level. 

In our numerics, the tight-binding Hamiltonian matrix 
is used to obtain the retarded Green's function of the 
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scattering region via the equation 

(E + irj -n-?:)G R = 1, 



(1) 



where the coupling to the leads is described by a self- 
energy S. It is then a straightforward task to obtain the 
scattering matrix from the Fisher-Lee formula [8j]: 



s = -I + 2 V ^ r * t P t G i? P* v / v- 



(2) 



Here v is a diagonal matrix containing the propagation 
velocities of the modes, ^ mn — \J w+ 1 sm w+i > wnere m 
numbers the sites in the leads and n is the mode index, 
and P mn is a matrix that contains the hopping ampli- 
tudes from graphene site m to lead site n. 

Once the scattering matrix has been found, we use it 
to calculate the linear conductance [t| 
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and the zero-frequency cross correlators 



S af 3 = Gq Tr [l^V] 



|/x«5-/x 7 |. (4) 



Here Gq = 2e 2 /h and ^ are the potentials in the leads. 
For noninteracting fcrmions, cross-correlations between 
different leads, a ^ 0, are always negative Q. Note 
that we defined the conductance such that it corresponds 
to linear response of the current I a in lead a when the 
potential (ip — —eVp in lead [3 is slightly varied and the 
potentials of all other leads are kept constant. 

According to Eq. ((4|, the cross-correlators are depen- 
dent on both direct processes coupling only the two leads 
(when 7 and 8 are both either a or /?) where the corre- 
lators are measured, or indirect processes involving also 
other leads. These two types of contributions can be sep- 
arated in an exchange Hanbury Brown- Twiss experiment 
d, • There the noise generated by one current source 
alone is compared to the noise generated by two current 
sources. Noise correlations between leads 1 and 2, de- 
fined in the following by S — ~S\2, are studied in three 
different cases: In experiment A voltage V is applied to 
reservoir 3, whereas in experiment B voltage V is applied 
to reservoir 4. In experiment C voltage V is applied to 
both reservoirs. The biasing scheme is shown in the in- 
set of Fig. El Classically, Sc — Sa + Sb, but quantum 
mechanical interference effects give rise to an exchange 
term, AS = S c - S A - S B - The latter is given by 0] 



AS = 
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In general, the exchange term can have either sign. Neg- 
ative AS means that exchange effects suppress noise 
whereas positive AS implies enhanced noise. Moreover, 
a finite AS* can be present only if there are modes that 



couple simultaneously all four leads so that the corre- 
sponding elements for all the scattering matrices above 
are non-zero. As we show below, the latter is the case in 
our system for evanescent modes close to the Dirac point, 
but not for the propagating modes far away from it 11|. 

We apply the multiprobe scattering formalism to a 
cross-geometry formed of a graphene sheet at the inter- 
section of four metallic leads. In our practical implemen- 
tation we use a graphene sheet with some 8500 lattice 
points (unless specified otherwise). The square lattice 
leads are matched to the graphene sheet so that the lat- 
tice constant in the leads connected to a zigzag edge is 
<2l = a and in the leads connected to an armchair edge 
ol = a/\/%, where a 2.46 A is the lattice constant 
in graphene. The Fermi level in the leads is chosen to 
correspond to half-filling of the band. This allows for a 
description of good contacts Q. 

Numerical simulations such as ours are often prone 
to effects related to the finite size of the simulated lat- 
tice. Such finite-size effects are relevant in the study 
of nanoribbons [l2| , but in typical experimental samples 
with dimensions in excess of 100 nm these effects are 
washed out. In our simulations, these finite-size effects 
are related to the change of the number of propagating 
modes inside graphene, and therefore lead to rapid os- 
cillations of the calculated quantities vs. gate potential. 
Moreover, in g raphene with zigzag edges there forms an 
edge state [13j which behaves differently from the rest of 
the states. In most experimental systems the effect of this 
state is fairly small, but it affects some of our numerical 
results. In our numerics, two of the edges are always of 
the zigzag type and two of them are of the armchair type. 
To take into account the effect of the type of the edges at 
the graphene-lead interface, we consider two alternative 
geometries, depicted in Figs. Ufa) and (b). 

The cross-conductances for a square sheet of graphene 
as a function of the chemical potential are shown in 
Fig. [21 The direct conductances obey the approximate 
symmetry Gi2(+A*g) ~ G^— /ig) and are quite small at 
the Dirac point /iq = 0. The "skew" conductance G13 is 
much larger at the Dirac point. It increases with increas- 
ing |^g| as well, but with a slope that is roughly half of 
that in the direct conductance. The evanescent modes 
are highly localized at the graphene-lead interfaces and 
therefore the coupling between adjacent leads is strong 
at the Dirac point, leading to a relatively high G13. 

We have also studied the geometry dependence of the 
cross-conductances by varying the size and aspect ra- 
tio of the graphene sheet. Figures Elja) and (b) show 
the cross-conductances as a function of the size of the 
graphene sheet. The aspect ratio was kept as close to 
unity as possible. For the evanescent modes (a) the di- 
rect conductances are constant, G w 0.15Gq, as can be 
expected from the pseudodiffusive model. The skew con- 
ductance increases with size, and approaches asymptot- 
ically the value G13 = 2Gq, one conductance quantum 
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FIG. 2: (color online) Cross-conductances in setup (a) of 
Fig. [T] as a function of chemical potential /j,g in graphene. 
Dashed lines show the predictions from the ray-optics pic- 
ture, Eq. (O, shifted by a constant to take into account the 
contribution of the evanescent modes. Note that the slope 
dG/dfic of these predictions is roughly in accord with the 
numerical results. 
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FIG. 3: (color online) Cross-conductances as a function of the 
size of the graphene sheet. In (a) and (b), the aspect ratio 
W/L = 1 is kept constant while the number of lattice points 
is increased, (a) is plotted at the Dirac point (evanescent 
modes) and (b) for /xgA ~ 0.5 (propagating modes). In (c) 
and (d) we plot the cross-conductances as a function of the as- 
pect ratio at the Dirac point and for na/t = 0.5, respectively. 
In (c) the graphene sheet had a fixed size of WL = 3600a 2 
and in (d) we chose W = 60a. In (d) the dashed lines show 
the prediction from Eq. (j6}. 



for graphene with spin and valley degeneracies. The de- 
viation of G13 from 2Gq is thus a finite-size effect, pre- 
sumabiycaused by the quasibound state near the zigzag 



edges [13} . For the propagating modes (b) all conduc- 
tances increase linearly with size, in line with the linear 
increase in the number of modes. 

Fig.[3Jc) shows the cross-conductances as a function of 
the aspect ratio of the graphene sheet for /j,q jt — 0. The 



size of the graphene sheet was kept constant. The direct 
conductances scale as G12 ex W/L and G34 oc L/W as 
expected from Ref. [l| with a conductivity of Ae 2 /nh. The 
skew conductance is quite insensitive to the aspect ratio, 
deviating from its bulk value only at very small or large 
aspect ratios where finite-size effects play a greater role. 

For large /iQ we can assume that all modes are prop- 
agating. In this case we can estimate the relative mag- 
nitudes of the direct and skew conductances with the 
ray optics model: In the ballistic limit all modes are 
ideally transmitted, and the terminal from which the 
electron exits is determined by the angle of propaga- 
tion, kpsinO = k y , where hp — 2^iq/ {\/Zta). For large 
structures we can convert the summation over transverse 
modes in the Landauer formula to an integral and obtain 
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Figure [Hd) shows the cross-conductances as a function 
of the length of the graphene sheet for /icA = 0-5. The 
sum of all cross-conductances is constant for large sam- 
ples and agrees with Eq. ((6|), confirming that all prop- 
agating modes are ideally transmitting. This also im- 
plies that evanescent modes, not taken into account in 
Eq. ©, do not contribute to conductance. The individ- 
ual cross-conductances deviate somewhat from the ray- 
tracing model, however, but the slope seems to be correct 
for large samples. This deviation is probably due to the 
relatively small size of our structure, where the replace- 
ment of the sum by an integral is not fully justified. 

We now turn to the noise correlations in the four probe 
setup. We first assume that current is driven between ter- 
minals 1 and 2 and the potential in terminals 3 and 4 is 
kept floating so that no average current flows in them. 
Due to symmetry the floating potential lies halfway be- 
tween the potentials in terminals 1 and 2. 

We define the local Fano factor with Fi oca i = Sn/2el, 
where Sn is the noise autocorrelator (shot noise) in ter- 
minal 1 and / the average current flowing from 1 to 
2. Similarly, the nonlocal Fano factor is defined with 
-Fnoniocai = S^/2eL. Symmetry again dictates that 
S\x = S22 and 5*33 = 5*44. These Fano factors are shown 
in Fig. [4] for both setups (a) and (b) as a function of the 
chemical potential in graphene with W = L. Each of 
them exhibit a peak of F « 0.37 at the Dirac point and 
decrease with increased doping of the graphene sheet. 
The oscillations at finite chemical potentials signal the 
appearance of propagating modes in the graphene sheet 
and their period is of the order of level spacing. The 
nonlocal Fano factors decrease more rapidly and exhibit 
weaker oscillations. Shot noise also obeys roughly the 
symmetry F^(^q) = F^(— (j>g), similar to the cross- 
conductances. Upon increasing the value of W/L the 
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FIG. 4: (color online) Local and non-local Fano factors as a 
function of the chemical potential in graphene. Lines with 
(a) or (b) in the legend refer to the corresponding setups in 
Fig. [T] Inset shows the dependence of the Fano factors vs. 
aspect ratio W/L at the Dirac point. 



local Fano factor of undoped graphene stays unchanged 
in both setups, but the nonlocal Fano factor decreases, 
reaching 0.1 at W/L = 10 (see inset of Fig. [J}. The de- 
viation of -Fiocai from 1/3, predicted in [l[ and observed 
in 



14j for two-terminal samples, is probably due to the 
effective inelastic scattering caused by the presence of the 
additional terminals 3 and 4. 

Figures E)Ja) and [HJb) show the exchange Hanbury 
Brown-Twiss cross correlations for setups (a) and (b), 
respectively. Near the Dirac point exchange correction is 
negative, and of similar magnitude with the classical con- 
tribution. For this reason two-terminal noise is almost 
identical to one-terminal noise, Si2/e|V|C?Q « —0.1. 
This finding is similar to the case of disordered box (lB | 
where transport is diffusive. In ballistic graphene, this 
result is due to the fact that evanescent states couple 
all the terminals. For larger no the exchange correction 
is almost vanishing (up to a one-channel interference ef- 
fect). The only differences between the two geometries 
are the magnitudes of the noise cross-correlations at in- 
creasing chemical potentials, which in setup (b) grow 
roughly quadratically versus the linear increase in setup 
(a). The negative exchange correction near the Dirac 
point seems to be quite robust, appearing also with non- 
square sheets and different interface transparencies. 

In conclusion, we have calculated cross-conductances, 
Fano factors and noise cross-correlations in graphene by 
applying a numerical tight-binding model. With these 
results we have pointed out that the semiclassical ray- 
optics picture usually valid in large ballistic normal con- 
ductors is also valid in doped graphene, but cannot be 
used to describe the behavior of the evanescent modes. 
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FIG. 5: (color online) Noise cross-correlation S = — S12 be- 
tween leads 1 and 2 as a function of chemical potential in 
graphene /j,g- (a) and (b) correspond respectively to setups 
(a) and (b) in Fig. [T] Inset shows the biasing scheme for the 
HBT experiment. 
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